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ABSTRACT 


A theoretical  method  is  derived  for  predicting  trim  angle  and  speed  coeffi- 
cient at  the  inception  of  porpoising  of  prismatic  planing  hulls.  Although 
equations  are  derived  for  the  surge,  pitch,  and  heave  degrees  of  freedom,  it  ir 
seen  that  the  effect  of  surge  is  small  at  ordinary  operating  trim  angles. 

Comparisons  of  theoretical  predictions  with  existing  experimental  data  on 
coupled  pitch  and  heave  porpoising  show  reasonably  good  agreement  for  a 
wide  range  of  speed  coefficients,  load  coefficients,  and  deadrise  angles.  The 
theoiy  may  also  be  used  for  estimating  the  natural  frequencies  and  damping 
characteristics  of  prismatic  hulls  in  the  stable,  high-speed  planing  range. 

ADMINISTRATIVE  INFORMATION 

This  investigation  was  authorized  and  funded  by  the  Naval  Sea  Systems  Command 
(SEA  035)  under  the  General  Hydrodynamics  Research  Program,  SR-023-0101,  Work  Unit 
1-1562-002. 


INTRODUCTION 

Porpoising  is  an  instability  in  pitch  and  heave  experienced  by  planing  craft  traveling  at 
high  speeds  on  calm  water.  It  has  been  known  to  lead  to  such  violent  motions  as  to  cause 
many  serious  boating  accidents.  With  constantly  increasing  boat  speeds,  this  phenomenon  is 
becoming  more  and  more  of  a problem  to  planing-  boat  designers. 

Perhaps  the  first  attempt  at  treating  this  problem  analytically  wis  made  by  Perring,1  who 
developed  a theory  for  porpoising  based  on  low-aspect-ratio  wing  theory.  The  practical 
application  of  this  theory  was  unsuccessful  since  the  theory  was  oversimplified.  Since  then,  a 
great  deal  of  experimental  work  has  been  done  concerning  porpoising  for  water-based  aircraft 
and  planing  boats,  and  more  recently  additional  attempts2,3  at  developing  a theoretical  treat- 
ment have  met  with  varying  degrees  of  success.  Perhaps  the  only  systematic  experimental 
investigation  for  planing  boats  was  done  by  Day  and  Haag4  on  prismatic  bodies.  These  bodies 
comprised  a wide  range  of  design  parameters.  The  results  of  this  work  have  been  widely  used  as 


'perring,  W.G.A.  and  H.  Glauert,  “Stability  on  lhe  Water  of  a Seaplane  in  the  Planing  Condition,"  Aeronautical  Research 
Council.  TR  Vol.  42  (Sep  1933).  A compiele  ! uling  of  reference*  is  given  on  pages  60  and  61. 

2lutowski,  R.N.."A  Computer  Program  for  Various  Performance  Aspects  of  Planing  Craft,"  Thesis  submitted  to  Stevens  Institute 
of  Technology,  Castle  Point,  Hoboken,  N.J.  (1973). 

^Paync,  P R..  “Coupled  Pitch  and  Heave  Porpoising  Instability  in  Hydrodynamic  Planing,"  Journal  of  llydronautics,  Vol.  8, 

No.  2 (Apr  1974). 

*Day,  J P and  R J.  Haig,  "Planing  Boat  Porpoising"  Thesis  Submitted  to  W;bb  InslHute  of  Naval  Architecture,  Gien  Cove, 

Long  Island,  N.Y.  (May  1952). 
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a guide  in  estimating  the  porpoising  limits  of  planing  hulls.5  Although  this  is  a reasonable 
empirical  approach,  it  does  not  contribute  much  to  a good  theoretical  understanding  of  the 
problem.  Such  an  understanding  is  required  to  determine  the  effects  of  variations  in  hull 
parameters  of  practical  boatr  as  well  as  to  evaluate  innovative  ideas  for  prevention  of 
porpoising.  Furthermore,  it  is  important  for  providing  a tool  for  estimating  the  effects  of 
design  parameters  on  natural  oscillation  frequencies  and  damping  characteristics  of  the  boat, 
since  these  characteristics  play  a dominant  role  in  dynamic  behavior  in  a seaway.  In  fact,  the 
theoretical  approach  derived  herein  has  been  used  ;n  development  of  a linearized  theory  for 
predicting  the  motions  of  planing  boats  of  arbitrary  deadrise  angle  in  waves.^ 

Although  the  theory  has  been  developed  for  prismatic  planing  hulls,  it  appears  to  be 
suitable  as  a guide  for  predicting  effects  of  parameter  variation  on  porpoising  of  practical 
planing-boat  configurations.  Furthermore,  the  methods  used  lend  themselves  to  direct 
extension  to  a theory  for  nonprismatic  hulls. 

STABILITY  EQUATIONS 

Stability  equations  for  the  longitudinal  motions  in  surge,  pitch,  and  heave  are  derived  in 
Appendix  A.  To  make  comparisons  with  the  large  quantity  of  existing  porpoising  data  from 
towed  models,  provision  was  made  to  include  the  effect  of  the  tow  force.  Since  the  models 
were  towed  at  constant  speed,  the  surge  equation  played  no  role  in  determination  of  the  boat 
model  stability;  only  the  coupled  pitch  and  heave  equations  were  needed  to  investigate  the 
problem.  However,  it  is  believed  that  the  results  obtained  in  this  manner  are  generally 
applicable  to  boats  with  all  three  degrees  of  freedom,  since,  as  shown  in  Appendix  A,  the 
magnitude  of  the  stability  derivatives  in  the  surge  equation  are  considerably  smaller  than  those 
in  the  pitch  and  heave  equations  for  most  cases  of  interest.  We  therefore  used  the  following 
nondimensional*  linearized  stability  equations  for  the  dynamic  heave  force  and  pitch  moment 
equilibrium,  respectively.** 


■’SjvitAv,  [>,,  "Uydrody ramie  Design  of  Planing  Hulls,"  Marine  Technology  (Oct  1964). 

^Martin.  M.,  "Theoretical  Determination  of  Motion  of  High-Speed  Planing  Craft  in  Wavet,"  DTNSRDC  Report  76-0069 
(Apr  1976). 

•The  prime  symbol,  normally  used  to  denote  a nondimenrional  quality,  is  omitted  for  convenience. 

''•The  effect  of  the  surge  degree  of  freedom  is  readily  determined  from  the  stability  equations  for  surge,  pitch,  and  heave  as 
derived  in  Appendix  A. 
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(Z"  - m)z  + Z^z  + Zzz  + Zq6  + Zg0  + Zg0  = 0 (1 ) 

M^z  + M^z  + Mzz  + (Mg  - ly)0  + Mg0  + Mg0  “0  (2) 

These  equations  describe  the  motion  relative  to  fixed  horizontal  and  vertical  axes,  0 x and  0 z, 
along  and  at  right  angles  to  the  direction  of  motion;  see  Figure  1.  The  origin  0 was  taken  at 
the  boat  center  of  gravity  and  moves  with  the  constant  reference  speed  U of  the  boat.  The 
symbols  Z and  z represent  the  verucal  force  and  displacement,  respectively,  at  0,  positive 
in  the  down  direction.  The  symbols  M and  6 are  the  pitch  moment  and  angular  displacement 
perturbation,  respectively,  with  respect  to  the  origin  0,  positive  in  the  sense  of  bow  up.  The 
coefficients  of  the  variables  z,  z,  z,  0,  0,  and  0 are  the  stability  derivatives;  e.g.,  Zg  is  the 
nondimensional  linearized  rate  of  increase  of  vertical  force  Z with  nondimensional  angular 
acceleration  0. 

The  stability  derivatives  have  been  derived  for  constant  deadrise  planing  hulls  in  Appendix 
A ;n  terms  of  the  geometric  and  operational  characteristics  of  the  boat.  They  were  derived  on 
the  assumption  that  the  craft  could  be  treated  as  a slender  body  with  an  empirical  three- 
dimensional  correction.  Because  of  the  high  Froude  number  range  of  operation  and  the  low 
aspect  ratio,  wavemaking  and  unsteady  lift  effects  were  assumed  negligible. 

The  velocity  and  acceleration  derivatives  (Appendix  A)  are 

Z £ - cos2  tJ" h'  ds'  (3) 

Zi  = -2  <^(X)  jis'  cos3  t (4) 

Zg  = - s')  ds'  (5) 

Zg  =-2^(X)fis'Xg  cos2  t (6) 

M”  = <p(X)  cos  rj  ji'(a'  - s' )*"  ds'  (7) 

M-  = 2^(X)cos2  v fa' (8) 

= -s')2  ds'  (9) 

Mg  = -2*(X)  cos  Xg2  + J \ u'(a'  - s')  ds'^  (10) 

where  r = equilibrium  trim  angle 

X = mean  wetted  length-to-beam  ratio 

<J>(X)  = three-dimensional  correction  factor 

li  = nondimensional  sectional  added  mass 

t i ' = nondimensional  sectional  added  mass  at  transom 
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\g  = nondiinensional  distance  from  transom  to  center  of  gravity 
s'  = nondimensional  distance  from  foremost  wetted  point  on  keel 
to  any  boat  section;  see  Figure  20 
a'  = value  of  s'  at  boat  center  of  gravity 


The  integrations  have  been  taken  over  the  wetted  length  of  the  boat.  Expressions  for  the 
sectional  added  mass  distribution  are  derived  in  Appendix  A. 

The  so-called  static  derivatives  are  obtained  directly  from  the  expressions  for  the  steady 
state  Zs'  force  and  Mg'  moment.  Since  both  are  functions  of  X and  r,  and  X = X(z,  r),  we 
have 


Zz  = 

Mz  = 
M0  = 


Mg  3X 

ax  bz 

(ID 

*J£  Ms  ax 

br  + aX  dr 

(12) 

™ s ax 
ax  3z' 

(13) 

3Ms'  3Mg'  ax 
br  3X  3r 

(14) 

The  steady  state  Zg'  force  which  is  the  negative  of  the  lift  force  (Appendix  A)  is 


Zc'  = ' sin  r cos^  r + ^ X^  sin  t cos  t -X  Cf  sin  t/ cos  0 

3 5 / 1 


(15) 


where  P = deadrisc  angle 

Cy  = U/\/"gF  = speed  coefficient 
Cf  = skin  friction  coefficient 


The  first  term  is  the  dynamic  lift  on  the  hull;  the  second  term  is  the  hydrostatic  lift,  and  the 
last  term  is  the  vertical  component  of  the  skin  friction  force,  assumed  to  act  parallel  to  the 
keel  line. 

The  expression  for  the  steady  state  moment  Ms'  about  the  boat  center  of  gravity 
(Appendix  A)  is 


■ zh  x (x«  - Hr)'  TS  f(ki • 'I ) + FWS  «2> 


Mg'  - <p(X)  sin  r cos  r 
C 


(16) 


4 


ft 

£ 


where 


Tg  = nondimensional  towing  force  = drag 
FWS  = nondimensional  aerodynamic  drag 
f(kj,ej)  = moment  arm  of  towing  force 
f(l<2,  e2>  = moment  arm  of  aerodynamic  force 

The  equation  for  Tg  (Equation  (90)  of  Appendix  A)  is 


tan  t + 


XCf 

cos  0 cos  r WS 


(17) 


The  first  term  in  Equation  (16)  is  the  hydrodynamic  moment,  the  second  term  is  the  hydro- 
static moment  and  the  third  term  is  the  moment  due  to  skin  friction.  The  last  two  terms  are 
the  moment  due  to  the  tow  force  and  the  aerodynamic  drag. 

The  steady  state  values  of  X and  r used  in  Equations  (3)  through  (14)  are  determined 
from  the  following  equations  of  planing  equilibrium. 

W'  + Zs'  = 0 (18) 

Mg'  = 0 (19) 

These  equations  are  solved  by  an  iteration  process  described  in  Appendix  A. 

The  solution  to  Equations  (I)  and  (2)  are 


a.'t' 


_ f.t 
O^y  t 


z = Zj  e 1 + z2  e * + 


a*  Y o?Y 
6=6  j e 1 + &2  e 1 + 


(20) 

(21) 


where  z. j . z-^  * * • 0],  #2’  * * ’ are  constants  which  depend  on  the  initial  conditions.  The  o' 
terms  determine  the  character  of  the  time  history  response  of  the  boat  to  any  small  disturbance. 
Four  values  of  o'  are  obtained  from  the  roots  of  the  resulting  characteristic  equation. 


A o'4  + Bo'3  + Co'2  + Do'  + E = 0 


(22) 


where  A = (Z-z-  - m)  (Mg  - Iy ) - MgZ$ 

B = Zg(M0  - Iy)  + (Z"  - m)M^  - M^Zg  - M^Z^j 

C = Zz(M y - ly)  + Zaftig  + (Zj  - m)M0  - M»Z0  - M^Z#  - MzZ$  (23) 
D = ZzMg  + ZgMg  - M;Zg  - MzZ^ 

E =ZzM0-MzZ0 

The  roots  of  these  equations  may  be  real  or  complex  conjugate  pairs.  In  either  case,  it  is 
seen  from  Equations  (20)  and  (21)  that  if  any  root  has  a positive  real  part,  the  transient 
response  increases  without  limit,  and  the  boat  is  considered  unstable  in  the  linear  sense. 
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In  general,  a complex  pair  of  roots  represents  an  oscillatory  mode,  e.g.,  for  the  root  pair 
o'  = oR'±ioj',  the  z'  response  is 

__  t .1 

z' = e ^ (C,  cos  o,'t'  + C,  sin  Oj't') 

where  Cj  and  C2  are  real  constants  which  are  determined  by  the  initial  conditions.  The 
magnitude  of  the  imaginary  part  of  the  root  Oj'  is  the  nondimension al  natural  frequency  of 
the  modal  motion.  In  dimensional  form  the  natural  frequency  and  period  are 


°I 


rad/sec 


(24) 


T 

T = — sec 
°\ 


(25) 


The  effect  of  the  real  part  of  the  root  oR'  may  be  illustrated  by  computing  the  time  for  a 
disturbance  to  either  halve  or  double  itself  in  magnitude.  Thus,  if  oR'  is  negative,  the 
envelope  of  the  disturbance  will  be  halved  when 


On  t Op t 

e k = e K = 


= 1/2 


It  follows  that  the  time  for  the  disturbance  motion  of  each  mode  to  halve  or  double  itself  is 


1 1 /2  or  t2  = 0.69 /oR  sec 


(26) 


Another  useful  measure  of  damping  of  oscillatory  modes  is  the  damping  ratio  £,  which  is 
directly  related  to  the  rate  of  decay  of  disturbance  oscillations. 

It  is  given  by 


* = - 


°R 

V °R  + °\ 


(27) 


In  the  vicinity  of  the  resonant  encounter  frequency  in  waves,  the  damping  ratio  is  also 
inversely  related  to  the  amplification  ratio  of  the  boat  response.  Values  of  £ between  0.6  and 
1.0  are  usually  considered  to  give  well-damped  modes.  Values  less  than  about  0.4  are  generally 
consideied  to  produce  underdamped  modes.  Although  the  forgoing  may  provide  a rough 
indication  of  the  vertical  plane  dynamic  characteristics  of  the  boat,  a dynamic  motions  analysis 
is  required  for  any  detailed  study. 
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d. 


X 


DETERMINATION  OF  PORPOISING  CONDITIONS 


Perhaps  the  only  systematic  experimental  investigation  of  porpoising  was  carried  out  in 
1952  by  Day  and  Haag.4  Measurement  by  other  investigators7,8,9  have  generally  been 
incidental  to  a broader  program  primarily  concerned  with  resistance.  The  experiments  of  Day 
and  Haag  were  carried  out  with  3.8-inch-beam,  prismatic  wood  models  towed  by  a light  line 
from  a point  slightly  forward  of  and  above  the  center  of  gravity.  Deadrise  angles  of  0,  1C.6, 
and  20.5  degrees  were  investigated.  The  load  coefficients  were  0.36,  0.48,  and  0.72.  For  each 
speed  the  boat  center  of  gravity  was  gradually  moved  aft  until  the  boat  porpoised.  The  trim 
at  which  this  occurred  defined  the  critical  trim  angle  and  provided  a point  on  the  trim  angle 
versus  speed  stability  boundary.  Values  of  the  mass  and  moment  of  inertia  for  practically  all 
the  models  at  porpoising  were  provided.  Thus,  it  was  possible  to  make  a theoretical  calculation 
for  each  test  condition.  Small  allowances  for  aerodynamic  effects  and  chine  comer  radius 
were  made.  These  are  described  in  Appendix  B. 

Calculations  of  each  of  the  four  stability  roots  for  each  test  condition  were  obtained  from 
Equation  (.22).  In  the  vicinity  of  porpoising,  the  least  stable  root  was  complex  in  each  case, 
indicating  that  the  response  to  a disturbance  would  always  be  oscillatory.  Figure  2 is  a typical 
plot  of  the  variation  of  the  real  part  of  the  least  stable  root  oRJ'  with  the  nondimensional 
longitudinal  distance  Xg  of  the  center  of  gravity  from  the  transom  for  three  of  the  test  speeds 
with  P = 10.6  degrees  and  = 0.48.  It  is  seen  that  the  stability  roots  become  negative 
(stable)  for  values  of  X less  than  about  0.20  and  greater  than  about  0.85  for  the  cases  shown. 

O 

We  note  that  in  the  stable  region,  corresponding  to  the  small  Xg  range,  the  equilibrium  trim 
angles  are  much  higher  than  in  the  stable  region  of  the  high  Xg  range.  Also,  the  magnitude  of 
the  stability  root  in  the  small  Xg  range  is  usually  quite  small,  so  that  the  damping  of  the  boat 
oscillations  may  usually  be  expected  to  be  poor.  Taking  an  example  from  the  figure,  we  find 
from  Equation  (26)  that  the  time  for  a disturbance  to  damp  to  half  amplitude  at  Cv  = 2.67 
would  be  more  than  15  seconds  for  a boat  with  a 15-fcot  beam,  and  oR1'  = -0.01.  This 
would  be  intolerable  in  the  presence  of  even  very  small  disturbances.  On  the  other  hand,  it  is 
seen  from  the  same  figure  that  moving  the  center  of  gravity  forward  of  the  upper  value  at 
which  the  root  becomes  neptive  results  in  increasingly  negative  values  of  the  stability  root  and. 
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1 ndsma,  G.,  “A  Systematic  Study  of  Rough-Water  Performance  of  Planing  Boils,”  Davidson  Laboratory,  Stevens  Institute  of 
Technology,  Hoboken,  N.J.,  Report  1275  (Nov  1969). 

fi 

Clement,  E.P.  and  D.  Blount,  “Resistance  Tests  of  Systematic  Series  of  Planing  Hull  F orms,"  Transactions  Society  of  Naval 
Architects  and  Marine  Engineers,  VoL  71,  pp.  491  -561  (1963). 

4 Davidson,  K.S.M.  and  A.  Suarez,  “Tests  of  Twenty  Related  Models  of  V-Bottom  Motor  Boats  - EMB  Series  50,"  David  Taylor 
Model  Basin  Report  R-47  (1949). 
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therefore,  considerably  higher  damping  for  this  mode.  For  example,  for  the  same  speed  and 
Xg  = 1.5,  we  find  oR)'  = 0.3,  and  the  time  for  a disturbance  to  damp  to  half  amplitude  is  of 
the  order  of  only  one-half  second.  Since  the  calculated  value  of  o{j'  at  this  condition  is  about 
0.5,  the  damping  ratio  from  Equation  (27)  is  found  to  be  approximately  0.50;  whereas,  for  the 
previous  condition  it  was  only  about  0.02. 

An  alternative  plot  to  Figure  2 showing  the  variation  of  the  calculated  oR]'  with  trim 
angle,  for  a complete  range  of  test  speeds,  with  0 = 0 degrees  and  = 0.72,  is  shown  in 
Figure  3.  This  plot  is  typical  of  those  used  to  determine  the  critical  trim  angle  as  a function 
of  the  speed  coefficient.  Where  two  trim  angles  at  zero  root  crossing  are  shown,  the  smaller 
trim  angles  are  taken.  The  porpoising  boundaries  thus  determined  are  discussed  in  the  next 
section. 


COMPARISONS  WITH  EXPERIMENTS  ON  PRISMATIC  HULLS 

Figure  4 shows  the  variation,  with  the  speed  coefficient  Cy,  of  the  critical  trim  angle  of 
a prismatic  planing  configuration,  with  a deadrise  angle  of  zero  degrees,  for  three  values  of 
load  coefficient.  The  points  in  the  figure  represent  the  experimental  data.4  The  curves  were 
obtained  from  the  present  theory.  It  is  seen  that  the  theory  agrees  reasonably  well  with  the 
experiments  for  = 0.36  and  0.72.  For  the  intermediate  load  coefficient  = 0.48,  the 
theory  gives  intermediate  values  of  critical  trim  angle;  whereas,  the  data  for  this  condition  are 
close  to  those  obtained  at  = 0.36. 

Figure  5 is  a similar  set  of  graphs  for  the  prismatic  hull  with  a deadrise  angle  of  10.6 
degrees.  Here,  the  agreement  between  the  theory  and  the  experimental  results  over  the  whole 
speed  range  is  quite  good  for  = 0.48  and  0.72.  However,  for  = 0.36,  the  theory  tends 
to  underpredict  the  critical  trim  angle  about  1 degree  at  the  higher  speeds.  Part  of  this 
discrepancy  is  due  to  the  fact  that  the  chine  became  unwetted  at  the  transom,  and  the  theory 
does  not  account  for  this.*  The  discrepancy  tends  to  be  greatest  at  the  high-st  speeds  and  for 
the  lowest  loading  conditions.  The  region  where  the  whole  chine  is  out  of  the  water  is  indi- 
cated approximately  by  the  dashed  portion  of  the  curve. 

Figure  6 is  a similar  comparison  for  a hull  with  deadrise  angle  of  20.5  degrees.  Here  the 
theoretical  curve  tends  to  fall  slightly  lower  than  the  data,  at  the  lower  speeds,  and  underpre- 
dicts the  critical  trim  angle  by  about  one  degree  at  the  higher  speeds.  In  this  case  the  chine 
became  unwetted  sooner  than  for  the  boat  with  a 10.6  degree  deadrise  angle.  The  approxi- 
mate region  of  occurrence  is  indicated  by  the  dashed  portion  of  the  curves.  An  attempt  was 


•The  theory  is  being  modified  lo  include  these  effeett. 
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made  to  determine  the  effect  of  modifying  the  magnitude  of  the  lower  limit  of  integration 
in  part  of  the  added  mass  integral  as  discussed  in  Appendix  A.  This  had  the  effect  of  further 
improving  the  predictions  as  shown  by  the  broken  curves  in  Figure  6.  This  slight  alteration  to 
ihc  theory  had  practically  no  effect  on  the  predicted  porpoising  boundaries  for  the  hulls  having 
0-  and  10.6-degree  deadrise  angles. 

Day  and  Haag  found  that  by  plotting  the  measured  trim  angle  against  (Cyj/2)1^2  or 
v/T^/C y,  the  separate  curves  for  each  load  coefficient  collapsed  into  a narrow  band.  This 
result  was  also  found  from  the  theoretical  calculations.  Figures  7 through  9 show  comparisons 
between  the  experimental  data  and  the  calculations  plotted  in  this  manner.  Except  as  noted 
earlier,  it  is  seen  that  the  overall  agreement  between  the  theory  and  the  data  is  reasonably  good 
for  all  of  the  deadrise  conditions  investigated. 

For  each  critical  trim  angle  tc,  the  corresponding  critical  position  of  the  boat  center  of 
gravity  X , must  satisfy  the  steady  state  equilibrium  equations.  Thus  the  degree  of  agreement 
between  the  measured  and  predicted  magnitude  of  Xgc  is,  in  part  at  least,  a measure  of  the 
accuracy  of  the  steady  state  equations.  Figure  10  shows  a comparison  with  theory  of  the 
measured  variation  of  Xgc  with  speed  coefficient  for  each  of  the  conditions  investigated  by 
Day  and  Haag.  It  is  seen  that  the  best  overall  agreement  was  obtained  for  the  hull  with  zero 
deadrise  angle.  For  the  deadrise  angles  of  10.6  and  20.5  degrees,  the  calculations  were  about 
10  to  15  percent  lower  than  the  measurements. 

Tli is  result  seems  to  be  consistent  with  the  tendency  of  the  steady  state  Equations  (18) 
and  (19)  to  underestimate  the  steady  trim  angle,  thus  requiring  a more  aft  position  of  the 
center  of  gravity  to  obtain  a given  trim  angle.  Comparison  with  the  steady  state  trim  data  in 
Reference  7 of  a similar  formulation  recently  proposed  by  Brown10  exhibits  the  same  tend- 
ency to  approximately  the  same  degree.  However,  the  empirical  formulation  presented  by 
Savitsky5  does  not  show  this  tendency  and  gives,  on  the  average,  better  agreement  with 
measured  steady  state  trim  angles.  As  will  be  seen  later,  the  use  of  this  formulation  generally 
gives  larger  values  of  Xg  for  a given  value  of  r* 

Another  interesting  feature  of  Figure  10  is  that  the  effect  of  load  coefficient  appears  to 
be  small.  This  is  mainly  due  to  the  fact  that  the  more  highly  loaded  configurations  of  the 
Day  and  Haag  models  had  smaller  radii  of  gyration.  As  will  be  seen  later,  the  theory  shows 
that  these  two  effects  are  generally  compensatory  in  their  effect  on  the  critical  center  of 
gravity  position.  It  is  also  worth  mentioning  that  the  theoretical  values  of  the  location  of  the 

*(*l)rovrn,  P.W.,  “An  L xperimental  »nd  Theoretical  Study  of  Planing  Surfaces  with  Trim  Flaps."  Davidson  Laboratory,  Stevens 
Institute  of  Technology,  Hoboken,  N.J.,  Report  1463  (Apr  1971). 

•The  reasons  for  this  discrepancy  are  being  investigated. 
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hydrodynamic  center  of  pressure  and  center  of  gravity  were  within  5 percent  of  each  other  in 
nearly  all  cases.  As  one  might  expect,  this  was  the  case  throughout  the  stable  operating  rangt 
as  well. 

Figure  1 1 is  a plot  of  the  theoretical  nondimensional  mean  wetted  length  of  the  boat  at 
porpoising  against  speed  coefficient  for  the  various  conditions  investigated.  Unfortunately, 

Day  and  Haag  did  not  measure  the  wetted  length  of  their  models,  so  no  direct  comparison 
with  the  data  is  possible.  Average  values  of  observed  wetted  lengths  of  the  spray  sheet  at  the 
chine,  which  did  not  vary  much  with  speed,  are  shown  for  a rough  comparison.  As  one  might 
expect  the  spray  sheet  length  relative  tr  the  mean  wetted  length  becomes  smaller  with  increas- 
ing deadrise  angle.5 

For  a typical  configuration,  Figures  12  and  13  show  the  separate  effect  of  nondimensiona. 
loading  and  radius  of  gyration  ky/'b  on  the  magnitude  of  the  critical  trim  angle  tq  as 
determined  from  the  theory.  It  is  seen  that  rc  becomes  smaller  with  both  increasing  radius  of 
gyration  and  load  coefficient.  However,  the  effect  is  not  large  - especially  for  small  values  of 
CLb.*  it  will  be  recalled  that  the  higher  the  loading  the  smaller  was  the  radius  of  gyration  for 
the  models  investigated  by  Day  and  Haag.  This  effect  also  contributed,  to  some  extent,  to 
collapse  of  their  data  tor  different  load  conditions.  A similar  effect,  though  more  pronounced, 
may  be  shown  to  hold  with  respect  to  the  critical  nondimensioiL?l  value  of  the  longitudinal 
centei  of  gravity  position  and  mean  wetted  length.  Figures  14  through  16  show  that  the  value 
ot  Agc  and  Amc  increases  with  both  increasing  and  increasing  ky/b.  Here  these  effects 
account  almost  entirely  for  the  tendency  of  the  Xgc  versus  Cv  curves  (Figure  10)  and  the 
Amc  versus  Cy  curves  (Figure  1 I ) to  collapse  into  a narrow  band  for  the  various  loadings 
shown. 

Although,  the  calculated  values  of  Agc  tend  to  be  low.  Figures  14  and  15  nevertheless 
illustrate  some  interesting  trends.  It  is  seen  that  the  value  of  AgC  has  an  increasing  tendency 
to  reach  a limit  with  increasing  deadrise  angle  and  decreasing  load  coefficient.  Also  the  value  of 
AgC  becomes  smaller  with  decreasing  nondimensional  radius  of  gyration  and  increasing  deadrise 
angle.  It  thus  appears  that  high  deadrise  angle,  low  loading,  and  a small  radius  of  gyration  permit 
a more  aft  location  of  the  critical  center  of  gravity.  Furthermore,  it  appears  that  the  tendency 

Q 

for  the  curve  to  turn  down  at  the  higher  deadrise  angles  may  explain  the  observation  by  Stolz"  and 
others  that  “deadrise  surfaces  which  are  not  deep  into  the  porpoising  range  often  regain  stability 
at  higher  speeds.” 

Relatively  recently  Fridsma  reported  a rather  limited  amount  of  experimental  data  on 
the  porpoising  of  prismatic  hulls  with  deadrise  angles  of  10,  20,  and  30  degrees  and  with 

•Theory  and  experiment  ihow,  however,  that  natural  frequency  of  oscillatory  motion  decreaci  significantly  with  increasing 
moment  of  inertia. 
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various  load  coefficients.  Experiments  for  the  20  degree  deadrise  case  were  carried  out  with 
models  of  length-to-beam  ratio  of  4,  5,  and  6 so  that  different  values  of  ky/b  were  obtained 
for  the  same  in  some  cases.  Figures  (17a)  and  (17b)  show  a comparison  of  measured 
values  of  critical  trim  angle  made  by  Fridsma  with  those  computed  from  the  theory.  The 
theoretical  curves  are  basically  cross  plots  of  curves  such  as  those  in  Figure  12  for  various  load 
coefficients.  They  provide  a convenient  method  of  indicating  the  effect  of  variations  in  !<y/b, 
separately  from  Cy,  and  on  the  magnitude  of  tq.  It  may  be  seen  from  the  figure  that  the 
theory  predicts  rc  quite  well  in  most  cases  but  that  there  are  insufficient  data  to  verify  the 
variation  with  ky/b  predicted  by  the  theory. 

The  calculated  values  of  XgC  and  Xmc  in  these  cases  were  about  28  and  15  percent  lower, 
respectively,  than  the  measured  values.  The  former  is  more  than  twice  the  difference  computed 
for  the  Day  and  Haag  experiments.  As  noted  earlier,  this  was  felt  to  be  due  mainly  to  the 
fact  that  the  steady  state  equations  tend  to  underpredict  steady  state  trim  angle.  Values  of 
X werc  then  rccalculated  for  each  of  the  theoretical  values  of  rQ  from  the  following  steady 
state  equations,  adapted  from  Reference  5. 


X =o.75X — 

P (5.21  Cv2/X2)  + 2.39 


where  Xmc  is  defined  by  the  following  equations 

cu  = TC'-'(0.0I20  * 0.0055  Xmc5'2/Ov2) 

CLb  ’ rLo  - 0 “<’5  0 C u°  “ 


(28) 


and  Xp  is  the  nondimensional  distance  of  the  center  of  pressure  forward  of  the  transom.  The 
assumption  that  the  center  of  pressure  and  center  of  gravity  are  coincident  was  found  from 
the  theory  to  be  valid  to  within  a few  percent  in  nearly  all  eases.  The  use  of  the  above 
equations  resulted  in  better  agreement  with  the  experimental  values  when  the  theoretical  values 
of  rf  were  substituted.  Table  1 shows  a comparison  of  Xgc  thus  obtained  with  the  measured 
values.  Also  shown  are  comparisons  with  the  measured  mean  wetted  lengths.  It  is  seen  that 
the  calculated  values  are  within  a few  percent  of  the  measurements  in  nearly  all  cases. 


COMPARISONS  WITH  EXPERIMENTS  ON  NONPRISMATJC  HULLS 

Although  the  theory  and  data  presented  thus  far  are  for  prismatic  hull  forms,  they  show 
trends  which  are  close  to  those  observed  on  a variety  of  nonprismatic  models.  Two  principle 
sets  of  data8,9  are  available  for  comparison.  The  DTMB  Series  62  models  of  Clement  and 
Blount8  are  pure  monohedran  hulls  with  a transom-to-beam  ratio  between  0.64  and  0.80  and 
constant  aft  deadrise  angle  of  12.5  degrees.  Those  of  Davidson  and  Suare/9  have  a low  chine 


TABLE  1 - COMPARISON  BETWEEN  CALCULATED  AND  MEASURED 
VALUES  OF  X«c  AND  Xmc,  USING  COMPUTED  CRITICAL  TRIM 
ANGLE  AND  METHOD  OF  REFERENCE  5 FOR  MEASURED 
PORPOISING  CONDITIONS  OF  REFERENCE  7 


Deadrise 

0 

deg 

Load 

Coefficient 

Ca 

Speed 

Coefficient 

CV 

XD_  Calculated 

X„,„  Calcula  ed 
me 

Xgc  Measured 

Xm„  Measured 
me 

10 

0.912 

3.83 

0.86 

0.94 

0.608 

3.00 

0.90 

0.99 

0.304 

3.33 

0.93 

1.06 

0.304 

2.00 

0.94 

1.00 

20 

0.608 

3.89 

1.07 

1.17 

0.608 

2.73 

1.03 

0.97 

0.304 

2.66 

1.00 

1.12 

0.304 

2.98 

1.02 

1.00 

30 

0.609 

2.73 

1.02 

0.95 

0.912 

3.85 

1.07 

1.00 

line  and  warped  bottom  with  a transom-to-beam  ratio  of  0.88  and  a mean  deadrise  angle 
between  4.0  and  7.3  degrees.  In  view  of  the  approximate  nature  of  the  comparison  it  was 
felt  that  existing  calculations  for  the  10.6-degree  deadrise  configurations  would  give  a suffi- 
ciently good  representation  of  the  trends  measured  in  the  previously  mentioned  data. 

The  porpoising  boundary  described  by  Clement  and  Blount  was  presented  as  a plot  of 
Ctb/Xgc  versus  the  volumetric  Froude  number  F^  where 

Cn,  = 2 CA/CV2 
Fv  * U,/-vT75  = CVI(CA)'« 

It  is  seen  that  the  critical  trim  angle  does  not  appear  anywhere,  and  the  porpoising  boundary 
is  expressed  mainly  in  terms  of  the  critical  position  of  the  center  of  gravity.  This  method  of 
plotting  was  an  apparent  attempt  to  collapse  the  measured  values  of  Xgc  for  the  wide  range  of 
load  coefficients  investigated.  Figure  18  shows  a comparison  of  the  theoretical  calculati  ns 
with  the  Series  62  data.  The  theoretical  curves  are  constructed  from  the  calculated  curves  of 
Figure  10  for  0 = 10.6  degrees  and  C^  = 9.36.  0.48,  and  0.72.  The  C^  values  are  in  the 
approximate  range  of  most  of  the  data.  The  values  of  Gy }/X  given  in  Reference  8 were 

nondimcnsionalizcd  with  respect  to  the  beam  at  the  center  of  gravity,  while  those  shown  in 
Figure  18  have  been  based  on  the  mean  of  the  maximum  beam  and  the  beam  at  the  transom 
in  order  to  provide  a more  realistic  comparison  with  the  constant  beam  case.  It  is  seen  that 
the  theoretical  curve  follows  the  same  trend  as  the  data  but  is  about  10  to  15  percent  higher. 
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This  is  mainly  a reflection  of  the  difference  between  predicted  and  measured  values  of  X 

6C 

shown  in  Figure  10.  It  may  be  seen  that  the  points  for  different  loading  conditions  will  fall 
on  separate  lines  and  that  the  trend  with  loading  for  each  model  is  close  to  that  given  by  the 
theory;  e.g.,  compare  Conditions  3,  4,  and  5 with  the  trends  shown  in  Figure  15. 

It  is  easy  to  shew  that  the  measured  effect  of  radius  of  gyration  exhibits  the  same  trends 
as  that  shown  by  the  theory;  see  Figure  14.  Although  the  radius  of  gyration  ky/b  of  the 
models  was  not  measured,  it  is  reasonable  to  assume  that  it  was  increasingly  larger  for  the 
models  with  larger  length-to-mean-beam  ratio  Lp/BpA.  By  comparing  the  trends  of  the  data 
for  different  models  with  approximately  the  same  load  coefficient  (Conditions  1 with  3 and  5 
with  6)  we  see  that  the  value  of  X..  follows  the  trends  with  kv/b  predicted  by  the  theory. 

Gement  and  Blount  found  that  the  slope  of  a straight  line  through  the  data  points  was 
about  -2.5.  This  implies  that,  for  a given  value  of  C^ 


since  C^  = 2 C^/Cy2.  From  Figure  10  it  is  seen  that  the  previously  described  equation  would 
fit  the  theory  and  experimental  data  for  the  prismai.c  hulls  quite  well.  Although  the  equation 
may  be  a good  approximation  for  the  data  shown,  it  is  seen  from  the  theoretical  curves  of 
Figures  14  and  15  for  ky  /b  constant  that  it  is  not  safe  to  assume  that  this  is  true  in  all  cases. 

The  data  of  Davidson  and  Suarez9  are  most  conveniently  compared  in  the  form  of  a 
porpoising  boundary  presented  by  Agnelli.11  This  is  shown  in  Figure  19  as  a plot  of  measured 
values  of  the  critical  trim  angle  rc  against  2/(Cyj /XgC).  Included  also  are  the  data  points  for  the 
Series  62  hulls.  These  were  recently  obtained  by  the  author  from  Mr.  Blount  and  are  not 
quite  the  same  as  those  shown  in  the  plot  of  Reference  1 1.  Shown  for  comparison  with  the 
trends  predicted  by  theory  are  straight  lines  drawn  through  points  obtained  from  the  calculated 
curves  of  Figures  5 and  10  for  the  case  of  0 = 10.6  degrees  and  C^  = 0.36,  0.48,  and  0.72. 
Although  many  of  the  points  do  not  fall  on  the  lines,  it  is  clear  that  the  trends  are  predicted. 
More  detailed  comparisons  with  the  individual  models  must  await  further  extension  of  the 
theory  to  include  more  general-type  hull  configurations. 

CONCLUSIONS 

A theoretical  method  has  been  derived  for  predicting  the  conditions  leading  to  porpoising 
in  the  surge,  pitch,  and  heave  degrees  of  freedom  of  prismatic  hulls  with  arbitrary  deadrise 
angle.  Comparisons  of  the  theory  with  the  porpoising  boundaries  measured  on  towed  models 

1 1 Agnelli,  J.C.,  "Evaluation  of  the  Trim  of  a Planing  Boil  al  Inception  of  Porpoising,”  presented  it  Spring  Meeting  cf  Society 
of  Nivil  Architects  and  Marine  Engineers,  Lake  Buena  Visla,  Fla.  (Apr  1973). 
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with  freedom  only  in  pitch  and  heave  showed  reasonably  good  agreement.  Since  no  porpoising 
data  with  all  three  degrees  of  freedom  are  available,  it  was  not  possible  to  check  the  theory  for 
this  case.  However,  from  a comparison  of  the  relative  magnitudes  of  the  coefficients  in  the 
surge  equation  it  appears  that  this  effect  is  small.  In  any  case  this  may  readily  be  investigated 
in  more  detail  with  the  three  degree  of  freedom  stability  equations  in  Appendix  A. 

The  stability  roots  obtained  from  the  characteristic  equation  provide  estimates  of  the 
dynamic  behaviour,  such  as  oscillation  natural  frequencies  and  damping  characteristics  of  the 
boat  in  the  stable  region. 

Although  the  theory  was  developed  for  prismatic  hull  forms,  it  appears  to  be  suitable  as 
a guide  in  estimating  the  porpoising  limits  and  dynamic  characteristics  of  more  conventional- 
type  planing  hulls,  as  well  as  the  effects  of  variations  in  several  of  the  parameters.  It  is  felt 
that  by  an  extension  of  the  analytical  methods  used  in  the  present  analysis  even  closer  agree- 
ment with  data  on  prismatic  hulls  and  conventional  boats  could  be  achieved,  and  a tool  suitable 
for  investigating  the  effects  of  detailed  design  modifications  could  be  obtained. 
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APPENDIX  A 

DERIVATION  OF  PORPOISING  STABILITY  EQUATIONS 

It  is  assumed  that  the  planing  craft  has  a prismatic  hull  of  constant  deadrise,  is  moving  at 
constant  speed  parallel  to  a calm  water  surface,  and  is  free  to  perform  small  perturbation 
motions  in  pitch,  heave,  and  surge  about  its  steady  equilibrium  attitude.  Since  the  theory  is 
concerned  mainly  with  the  high-speed,  low-aspect  ratio  condition,  it  is  assumed  that  the  craft 
may  be  treated  as  a slender  body  with  an  empirical  three-dimensional  correction,  and  unsteady 
effects  are  small.  The  kind  of  analysis  to  be  used  was  first  used  in  1924  by  Munk1^  and  later 
by  Jones13  in  connection  with  the  analysis  of  airships  and  slender  wings,  respectively.  More 
recently  this  method  has  been  generalized  by  Bryson1**  for  completely  submerged  slender- 
finned  missiles.  It  has  also  been  applied  to  the  problem  of  pure  translational  impact  of  sea- 
planes on  a calm  water  surface  by  Mayo13  and  others.16,1^ 

FORCES  DUE  TO  PERTURBATIONS  IN 
VELOCITY  AND  ACCELERATION 

The  flow  over  the  hull  is  assumed  to  occur  in  transverse  planes  which  are  fixed  in  space 
and  oriented  normal  to  the  keel;  see  Figure  20.  The  momentum  of  each  layer  of  water  trans- 
verse to  the  keel  is  jifds,  where  is  the  two-dimensional  added  mass  of  the  section  of  the  hull 
at  point  s,  interacting  with  the  section  of  the  flow  plane  of  length  ds,  and  f is  the  component 
of  the  velocity  of  the  body  normal  to  the  keel  at  that  point.  The  coordinate  s is  measured 
from  the  foremost  immersed  station  along  the  keel.  The  normal  force  on  the  section  ds  of 
the  hull  is  the  time  rate  of  change  of  the  momentum  of  the  layer  of  water  ds  at  s. 


i 

*Munk,  M.M.,  “The  Aerodynamic  Force*  on  Airship  Uulli.”  National  Advisory  Committee  for  Aeronautic!  Report  184  (1924) 

1 ^Jonei.  R.T.,  “Propertiei  of  Low-Aipect-Ratio  Win(ti  at  SpeeJs  Below  and  Above  the  Speed!  of  Sound,"  National  Advisory 
Committee  for  Aeronautics  Report  stj  ( j 946) . 

'^Bryson,  A.E.,  Jr.,  “Stability  Derivative  'r  a Slender  Missile  with  Application  to  a Wing-Body  Vertical  Trail  Configuration," 
Journal  of  Aeronautical  Sciencei,  Vol.  - No.  5,  pp.  29 / - 308  ( I9S !». 

1 *Mayo  W.L.,  “Analysis  and  Modification  of  Theory  for  Impact  of  Seaplanes  on  Water,"  National  Advisory  Committee  for 
Aeronautics  Report  8i0  (1945). 

1 ^Milwitaky,  B.,  “A  Generaliied  Theoretical  and  Experimental  Investigation  of  the  Motions  and  Hydrodynamic  Loads 
I xperienced  by  V-Bottom  Seaplanes  During  Step-Landing  Impacts,"  National  Advisory  Committee  for  Aeronautics 
TN  1516(1948). 

1 'Schnltrer,  E..  “Theory  and  Procedure  for  Determining  Loads  and  Motions  in  Chine-Immersed  Hydrodynamic  Impacts  of 
Piismatic  Bodies."  National  Advisory  Committee  for  Aeronautics  Report  i 152  (1953). 
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(29) 


iFd  = jW)4 

Both  /j  and  f will  in  general  be  functions  of  the  longitudinal  position  coordinate  x and  time  t. 
The  time  derivative  is  therefore 


d_ 

dt 


(30) 


where  U is  the  steady  state  speed. 

The  normal  hydrodynamic  force  over  the  entire  hull  is  obtained  by  integrating  Equatior. 
(29)  along  the  wetted  length  of  the  hull  lk  and  multiplying  by  a correction  factor  0(X)  to 
account  for  the  three-dimensionality  of  the  flow. 


F 


D 


d 

mj  ^oiftds 


(31) 


where  X is  the  mean  wetted  length  divided  by  the  beam.  A plot  of  0(X)  obtained  empirically 
by  Pabst18  is  shown  on  Figure  21.  The  integral  may  be  expressed  as  the  sum  of  a velocity 
term  and  an  acceleration  term. 


FD=*(X) 


r r8MM . . rcx  df . i 

U Jo  ^ J 


(32) 


The  longitudinal  and  heave  perturbation  velocities  and  accelerations  are,  respectively, 
denoted  by  u,  u,  z,  z.  The  pitch  angle  perturbations  are  0,  0.  From  Figures  1 and  20  we 
obtain  the  following  relationships: 


ds 

dx 

a £ = 

3x 


= - COS  T 


- sin  r 


— = u sin  r + z cos  t - 0 (a  - s) 


f = 0 


(33) 

(34) 

(35) 

(36) 


where  r is  the  equilibrium  trim  angle  of  the  boat,  and  a is  the  value  of  s at  the  tranverse  flow- 
plane  through  the  boat  ctnter  of  gravity.  From  these  equations  and  Equation  (30)  we  have  to 
the  first  order  in  the  perturbations 


**Pjbst,  W.,  “Landing  Impact  of  SeipUitei,”  National  Advisory  Committae  for  Aeronautics  TM  624  (1931). 
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f = U sin  r + u sin  r + z cos  r - 0 (a  - s)  (37) 

$■  = 2 U0  cos  r + u sin  r + z cos  r - 0 (a  - s)  (38) 

$-f§ HI- 

On  substituting,  these  equations  into  Equation  (32),  dropping  the  ‘•.eond  order  perturba- 
tion terms,  and  integrating,  we  obtain 


FD  " FDS  + A FD 


where 


FDS  = s‘n  r cos  r 

A Fjj=  (2Fds/U)(u  + z cot  7 + 0 Cg/sin  r) 

I • rfik  •• 

+ i^(X)v(ii  sir.  r + z cos  r)  / juds~0/  , 


fi(a  - s)ds 


The  term  /j$  is  the  value  of  the  sectional  added  mass  at  the  stem  and  1=1^-  a.  The  first 
term  FDS  is  the  steady  slate  hydrodynamic  normal  force,  while  the  remaining  terms  A FD  are 
the  linearized  force  contributions  from  surge,  pitch,  and  heave  perturbations  in  velocity  and 
acceleration.  The  factor  2 multiplying  the  z and  0 terms  in  Equation  (42)  arises  from  the  fact 
that  the  sectional  added  mass  is  not  only  a function  of  position  on  the  hull  but  also  of  depth. 
This  is  u presented  by  the  contribution  of  df/dt  of  Equation  (35)  to  d/a/dt  in  Equation  (39). 

The  hydrodynamic  moment  is  obtained  by  integrating  the  product  of  the  stripwise  force 
in  Equation  (29)  and  the  moment  arm  from  the  center  of  gravity  a - s. 

/*  E|£  J 

Mp  = <0(M  / (a  - s)  — (/af)ds  (43) 

-'o 

With  the  aid  of  Equations  (31  through  39)  we  f.nd  from  Equation  (43),  after  dropping  second 
order  terms  in  the  perturbations 

Md  = Mpjj  + AMq  (44) 

where 


Mno  = tf(X)  O sin  r cos  t 


M ds  - £g 


(u  + z.  cot  r) 


AMn  = 


2MDS 


U 


- v?(X)  2 U cos  r j^s  eg2  *r  /i  (a  — s)dsj  6 

rk 

+ v?(X)(u  sin  t + z cos  r)  / p(a  - s)ds 

~'o 

rk 

- >p{\)6 ] p(a  - s)2  ds 


(46) 


The  first  term  Mpg  is  the  steady  state  hydrodynamic  moment  about  the  center  of  gravity  of 
the  boat.  The  remaining  terms  AMp  are  the  linearized  contributions  of  the  surge,  pitch,  and 
heave  perturbations  in  velocity  and  acceleration. 

VELOCITY  AND  ACCELERATION  STABILITY  DERIVATIVES 

The  vertical  and  horizontal  components  of  the  force  stability  derivatives  with  respect  to 
the  velocity  and  acceleration  perturbations  are  the  coefficients  of  the  perturbation  terms  in 
the  vertical  and  horizontal  components  of  the  perturbation  force  AFq  in  Equation  (42). 

We  write  these  equations  in  nondimensional  form*  by  dividing  through  by  1/2  pU  b . Thus 
the  vertical  component  of  the  perturbation  force  equation  becomes  with  Z positive  down 


azd  = 


2 FRS'  cos  r ( u'  + z cos  r + O'  Xg/sin  r) 


- <p(A) 


(u  sin  r cos  r + z c 


z' cos ~ r)F  p' ds ' + O' f p'(a'  - s')ds'| 


(47) 


where 


Fqs’  = <p(A)  < sin  r cos  r 

\ - 

p'  = p/i  ph2 


O'  = 0 b /U  . etc. 


(48) 


•Nondimensiorul  quantities  are  represented  by  : prime  symbol. 
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Likewise  the  horizontal  component  with  X positive  forward  becomes 


AXq'  = AZq'  tan  t 


(49) 


The  corresponding  nondimensional  stability  derivatives  for  the  moment  equation  are  obtained 
in  an  analogous  manner  from  the  nondimensional  form  of  Equation  (46)  as 


where 


AMd'  = 2Mds'  (u'  + z'  cot  r) 


V?(X)  2 cos  r 


>\2*J  V<*' 

r\ 

+ i^(X)(u'  sin  r + z'  cos  r)J  n'(a'  - s')  ds' 

“'o 

- y V(a'  - s')2  ds' 


Mds'  = *<X>  f V(a'-s') 

''a 


2 ds' 


(50) 


(51) 


Typical  notation  for  the  nondimensional  stability  derivatives  are  shown  as  follows 
zd'  * ~ 2 FDS  Xg  cot  r 


/xk 

fx'  ds' 


*U  = Zu  tan  r 


M, 


)'  * - *(K)J  V(a'  - s')2  ds' 


Xi'  = Z i tan  r 


Mu'  = 2Mds' 
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SECTIONAL  ADDED  MASS  DISTRIBUTION 


The  next  step  is  to  find  an  expression  for  the  distribution  of  the  boat  sectional  added 
mass  /i  as  a function  of  s.  A general  theoretical  expression  for  this  quantity  is  not  available, 
even  for  prismatic  hulls  with  constant  deadrise.  However,  relatively  simple  approximate 
expressions  for  added  mass  have  been  used  successfully  in  the  past  in  the  analysis  of  hydrody- 
namic impact  of  prismatic  bodies.  For  the  sections  of  that  portion  of  the  body  with  the 
chines  above  the  water  surface,  the  sectional  added  mass  was  estimated  by  the  following 
equation 

M = ^~2~  W)2  (52) 

where 

f(0)  = ^-l  (53) 

(3  - deadrise  angle  in  radians 

This  expression  is  based  on  the  work  of  Wagner.19  The  quantity  ff(0)  is  the  radius  of  the 
semicircular  cylinder  representing  the  added  mass  of  the  section.  For  the  sections  of  the  hull 
with  the  chine  submerged,  the  following  expression  was  used. 

2 

M = (f(0)tan0)2  + B|b(f-fc)  (54) 

where  B is  a function  of  the  angle  of  deadrise,  and  b is  the  beat  beam.  The  first  term  is  the 
contribution  of  the  V-shaped  bottom  alone  at  the  instant  of  chine  immrsion.  This  is  obtained 
from  Equation  (52)  by  putting  f = fc  = b/2  tan  0.  The  second  term  is  an  estimate  of  the 
effect  of  chine  depth  as  suggested  by  Schnitzel12  and  is  based  on  the  theory  of  Bobyleff20 
for  infinite  immersion.  The  Bobyleff  function  B is  shown  in  Figure  22. 

The  expressions  for  m given  by  Equations  (52)  and  (54)  will  used  in  Equations  (40) 
and  (44)  to  determine  the  normal  force  Fp  and  moment  Mp.  However,  it  is  first  necessary  to 
define  the  range  of  keel  length  over  which  each  contribution  to  v is  valid.  It  is  clear  that  no 
single  location,  such  as  suggested  by  Equations  (52)  and  (54),  exists  at  which  the  effect  of 
chine  immersion  starts,  since  the  flow  is  much  more  complex  than  these  equations  suggest. 

This  is  especially  true  in  the  vicinity  of  chine  immersion.  However,  a practical  solution  to 
this  problem  has  been  made  possible  by  making  the  theory  for  the  steady  part  of  the  normal 
force  and  moment  consistent  with  the  large  amount  of  existing  steady  state  data. 

19 

Wagner,  H..  "The  Phenomena  of  Impact  and  Planing  on  Water.”  National  Advisory  Committee  for  Aeronautics  Translation 
1366,  ZAMM  Bd  12.  Heft  4,  pp.  193  - 215  (Aug  1932). 

^°Lamb,  H.  "Hydrodynamics,"  Sixth  Edition.  Cambridge  University  Press,  England  (1932). 
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By  assuming  that  the  planing  hull  normal  force  at  high  speeds  was  made  up  of  the  sum 
of  low-aspect-ratio  wing  lift  and  a cross  flow  drag  term,  Shuford21  arrived  at  the  basic  form 
of  an  expression  which  he  was  able  to  fit  very  well  to  data  obtained  by  many  investigators. 

These  data  cover  a range  of  trim  angles  between  2 and  30  degrees,  wetted  lengths  from 
one  to  seven  beams,  and  deadrise  angles  between  0 and  50  degrees.  His  expression  for  the 
normal  force  is  given  by 


where 


fds  =Tpu2s[2(Tf  "ai  sin  T cos  T (l  _sin 

+ Cq  c sin2  r cos2  r cos  (3 

Cqc  = cross  flow  drag  coefficient 

S = plant'orm  area  of  the  wetted  portion  of  the  hull  or  Xb2 
A = aspect  ratio  of  S,  i.e. 


(55) 


A = b2/S  = 1/X  (56) 

The  values  of  Cq  and  the  dependency  on  deadrise  angle  were  obtained  by  a fit  to  the  data. 
The  cross  flow  drag  coefficients  were  found  to  have  the  values  shown  in  Table  2. 21 

TABLE  2 - CROSS  FLOW  DRAG  COEFFICIENT21 

Section  shape CDc 

V-bottom,  constant  deadrise 1.33 

V-bottom,  horizontal  chine  flare 1.33  + 0.0147  0° 

V-bottom,  vertical  chine  strips 1.60  + 0.0147  (3° 


We  will  now  define  the  various  contributions  to  the  sectional  added  mass  along  the 
hull  by  the  following  equations,  which  are  somewhat  less  restrictive  than  Equations  (52)  and 
(54).  , 


~~  f(0)2  s2  tan2  r 

o < s < sc, 

(57) 

^ f(0)2  scI2  tan2  r 

Scl  < s < «k 

(58) 

T-  B b tan  r (s  - sc-i) 

sc2  < s < Ck 

(59) 

V 


^Shuford,  C.L.,  Jr.,  “A  Theoretical  and  Experimental  Study  of  Planing  Surface!  Including  Effects  of  Cross  Section  and  Plan 
I orm,"  National  Advtory  Committee  for  Aeronautic*  Report  1355,(1957). 
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where  we  have  made  the  substitution 


£ = s tan  r (60) 

The  sectional  added  mass  n at  any  section  is  simply  the  sum  of  the  contributions  at  that  section. 
The  principle  difference  from  Equations  (52)  and  (54)  is  that  scl  and  sc2  have  not  been 
assumed  to  be  known  in  advance.  Equations  (57)  and  (58)  are  the  contributions  from  the 
bottom  of  the  hull  to  the  chine.  The  quantity  scl  is  the  value  of  s at  the  point  where  the 
chine  is  effectively  immersed.  Equation  (59)  is  the  contribution  corresponding  to  the  secon  J 
term  in  Equation  (54),  and  sc2  is  the  value  of  s at  which  this  begins  to  grow. 

From  Equations  (57  through  59)  we  readily  find  that 

Ms  = *7  (f(0)  sc,  tan  r)2  Bb  tan  r (£k  - sc2)  (61) 

2 

J *Vds  = ~ (f(0)  scl  tan  r)2  (ik  -jsc|j  + ^ B b tan  r — — j--2--  (62) 


Substituting  Equation  (61)  into  Equation  (41)  gives 

U2 


Fds  = — r~  M0)ir(f(0)  scl  tail  tY  sin  r cos  r 


+ i^(X)  B b sin2  r (£k  - sc2)] 


(63) 


This  equation  becomes  identical  with  the  formulation  of  Shuford  in  Equation  (55)  provided 


= 


X 

1 + X 


1 

I + A 


(64) 


(f(0)  scl  tan  r)2  - ^ ( I - sin  0) 

(65) 

^(X)  B = CD  c cos2  t cos  0 

(66) 

ek  “ sc2  = X b 

(67) 

The  three-dimensional  correction  factor  of  Equation  (64)  has  a trend  similar  to  the  result 
obtained  by  Pabst  as  shown  by  Figure  21.  Also  the  variation  of  B with  deadrise  angle  is  seen 
to  be  proportional  to  cos  0 in  Figure  22. 

If  we  substitute  Equations  (64)  through  (67)  into  Equation  (63)  and  nondimensionalize, 
we  obtain 

F'ds’  = | ^ sin  r cos  r ( 1 - sin  0) 


+ C 


D.c 


X sin2  t cos2 


t cos  0 


(68) 
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Likewise  if  we  substitute  Equations  (61)  through  (67)  into  Equation  (45)  and  divide  by 
1/2  pU2  b3  we  obtain  the  nondimensional  dynamic  moment  about  the  center  of  gravity 

MDS  = Y7J  ^sin  r cos  r (1  - sin  0)  (xk  - 1 Xc,  - Xg) 

+ Cq  c sin2  r cos2  r cos  0 - Xg^  X (69) 

where  XcJ  = scl/b 

Xc2  = sc2^b 

xg  = yb 

Xk  = Ck/b 


DETERMINATION  OF  XcJ,  Xc2,  Xk 


Because  of  wave  rise  on  impact,  the  effective  depth  of  the  V-bottom  is  greater  than  the 
depth  relative  to  the  calm  water  free  surface.  Wagner1^  found  that  before  chine  immersion, 
the  effective  depth  was  greater  by  a factor  of  rr/2.  Therefore  the  effective  radius  c of  the 
semicircular  cylinder  of  water  representing  Ms)  is  taken  as 


_ ir  tan  r 
2 S tan  0 

Combining  with  Equation  (57)  gives 


(70) 


Ms) 


* ? 17,  20\2  tan2  /f|  f 

P2C  [('-T)  — J fors<sO 


(71) 


The  factor  in  brackets  represents  the  effect  of  deadrise  angle  on  the  sectional  added  mass. 
Shuford  found  that  this  did  not  correlate  well  with  data  for  angles  of  deadrise  greater  than 
25  degrees.  He  therefore  substituted  the  function  (1  - sin  0)  which  correlated  well  with  data 
to  0 = 50  degrees.  Making  this  substitution  in  Equation  (71)  and  substituting  into  Equation 
(65)  we  find  the  following  expression  for  c 

c = v/T^^  for  s = sc]  (72) 

This  equation  is  seen  to  be  independent  of  deadrise  angle.  On  substitut.,1-;  into  Equation  (70), 
we  obtain 


scl 


_ y/T  tan  0 


tan  r 


Xcl  b 


(73) 
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An  expression  for  the  mean  wetted  length-to-beam  ratio  X has  been  obtained  by  Brown1® 


from  extensive  photographic  observations  as 

X = 0.5  (Xk  + Xc)  + 0.03  (74) 

where  X£,  the  ratio  of  the  wetted  length  of  chine  to  beam,  is 

Xc  = Xk  - (0.57  + 0.001  0)(tan  0/(2  tan  r)  - 0.006  0)  (75) 

provided  Xc  > 1 The  last  term  in  Equation  (74)  is  an  allowance  for  stagnation  line  curvature. 
We  then  find  from  the  previous  equations  that 

Xc2  = 0.5  (0.57  + 0.001  0)(tan  01(2  tan  r)  - 0.006  0)  - 0.03  (76) 

since  according  to  Equation  (67) 

= X + Xc2  (77) 


Equations  (73),  (76),  and  (77)  completely  define  the  ranges  of  the  sectional  added  mass 
distributions  in  terms  of  0,  X,  and  t.  Both  X and  r are  obtained  from  the  steady  state  equili- 
brium conditions  to  be  discussed  in  a later  section. 

It  is  noted  that  the  value  of  Xc2  as  defined  by  Equations  (59)  and  (76)  is  smaller  than 
Xk  - Xc  which  corresponds  to  the  lower  limit  of  integration  at  which  f = f in  the  representa- 
tion of  the  second  term  of  the  added  mass  in  Equation  (54).  Although  it  appears  more 
reasonable  to  use  Xk  - Xc  for  the  lower  limit  of  integration,  the  best  fit  to  the  data1®  requires 
that  we  use  the  Xc2  of  Equation  (76).  In  the  numerical  analysis,  the  effect  of  using  Xk  - X£ 

in  place  of  Xc2  was  found  to  be  insignificant  except  for  the  20-degree-deadrise  case. 
Calculations  of  r£  using  Xk  - Xc  are  shown  by  the  broken  curves  in  Figure  6 where  it  is  seen 
to  produce  a small  improvement  in  the  fit  to  the  data. 

ADDED  MASS  FUNCTIONS 

The  added  mass  functions  used  n<  evaluating  the  stability  derivatives  may  now  be 
expressed  in  terms  of  the  hull  geometry  and  the  integration  limits  Xc],  Xc2,  and  Xk. 

The  nondimensional  sectional  added  mass  at  the  stern,  including  the  three-dimensional 
effects  is  readily  obtained  from  Equations  (61)  and  (64)  through  (67)  as 

#)(is'=2(P  + Q)  (78) 

where 

p = 7 TTx  ( 1 ~ sin  0’ 
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c 

Q = — y~  X sin  r cos  r cos  0 

The  nondimensional  added  mass  in  heave  is,  from  Equation  (62) 

*(A)/XV  ds'  = 2PX4  +QX 

■'o 

where 


(79) 


X4  - Xk  - jXc, 

The  nondimensional  first  moment  of  the  added  mass  with  respect  to  the  point  s = 0 is  readily 
obtained  with  the  aid  of  Equations  (57)  through  (59)  as 

•Xu 

...  I lit 

(80) 


*00 / V s'  ds'  = 2^Px2  + ^*3) 


where 


\ Acl 
x2  2 ~ 4 


*3 


\3  - >.23  A 2 - 1 2' 


3 \ ? / ac2 

The  nondimensional  second  moment  of  the  added  mass  with  respect  to  the  point  s = 0 is 

•Xu 


where 


SP(X)  J V s’2  ds'  = 2 ^Px5  + 


\!3  . Xk3  ' Xd3 
X5  = ~5  + 3 


(81) 


Xu4  -X 


*6 


c2 


-(^4 


lc2 


REMAINING  FORCE  AND  MOMENT  CONTRIBUTIONS 

To  completely  specify  the  forces  and  moments  on  a planing  craft,  it  is  necessary  to  add 
to  Equations  (40)  and  (44)  the  contributions  due  to  perturbations  in  pitch  and  heave  displace- 
ment, buoyancy,  skin  friction,  thrusters,  aerodynamic  factors,  and  towing  forces  and  moments. 
Although  the  thrusters  may  have  a significant  effect  on  the  equilibrium  trim  and  wetted  length, 
their  effect  on  the  stability  derivatives  is  probably  not  very  large.  In  any  event  these  effects 
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may  be  estimated  by  means  of  existing  techniques,22,22  and  will  not  be  considered  further 
here.  However,  the  effect  of  a towing  force  and  moment  will  be  included  in  order  to  facilitate 
comparison  with  towed  model  data. 

The  contributions  due  to  perturbations  in  pitch  and  heave  displacement  may  readily  be 
obtained  from  the  complete  expressions  for  the  steady  state  force  and  moment.  The  additional 
terms  required  to  complete  these  expressions  are  described  as  follows. 


BUOYANCY  FORCE  AND  MOMENT 


We  will  consider  only  the  cases  for  speed  coefficient  Cy  greater  than  0.5,  where  the  water 
breaks  clear  of  the  transom,  thus  fully  ventilating  the  backside  of  the  boat  to  the  atmosphere. 
For  this  case  the  hydrostatic  force  may  be  assumed  to  act  normal  to  the  keel.  The  following 
expression  was  found  to  fit  the  data  reasonably  well.  10,24 

Ffis'  = « X2  sin  t/Cv2  (82) 


where  k is  an  empirical  correction  factor  which  accounts  for  ventilation  effects  on  the  static 
pressure.  On  the  basis  of  preliminary  analysis  of  planing  boat  test  data,  a value  of  k of  0.7 
was  tentatively  suggested  by  Hsu.24  However,  from  recent  extensive  experiments  with  a 10- 
degree  deadrise  prismatic  planing  hull,  Brown10  obtained  the  best  agreement  with  the  data  by 
putting 

k = 0.624  (83) 


and  assuming  that  it  acts  at  one  third  of  the  mean  wetted  length  from  the  stem.  Since  the 
present  analysis  deals  with  prismatic  hulls  Equation  (83)  will  be  used  in  the  following.  The 
moment  about  the  center  of  gravity  is  clearly 

MBs'sFBs'(j-,'g)  <84> 


SKIN  FRICTION 

The  contribution  due  to  skin  friction  is  assumed  to  act  tangential  to  the  bottom  and  mid- 
way between  the  keel  and  chine  lines.  It  is  given  in  terms  of  the  mean  wetted  area  by 
Savitsky5 * * * * 


5 

“Midler,  J.B.,  “The  Prediction  of  Power  Performance  on  Planing  Cr»ft,”  Transactions  Society  of  Nival  Architect!  ind  Marine 

f nginee, »,  Vol.  74.  pp.  563  - 610  ( 1966). 

^Ribner,  H.S.,  "Propellers  in  Yaw,"  National  Advisory  Committee  for  Aeronautics  Report  820  (194V). 

^*Hsu,  C.C.,  "On  the  Motions  cf  High  Speed  Planing  Craft,"  Hydronautic*  Report  603-1  (May  1967). 
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(85) 


Fpg  = X Cf/cos  0 


where  the  friction  coefficient  Cf  is  given  by 


0.242 

V'q 


as  a function  of  Reynolds  number.  The  moment  about  the  center  of  gravity  is  clearly 


(86) 


TOWING  AND  AERODYNAMIC  FORCES 

Terms  similar  to  the  above  Equation  (86)  may  be  written  for  the  towing  and  aerodynamic 
force  and  moment  contributions.  These  are  assumed  here  to  act  in  a line  parallel  to  the  steady 
part  of  the  straight  line  motion.*  The  moments  about  the  center  of  gravity,  due  to  the  tow 
force  T<j'  and  wind  force  FWS«  respectively,  are 

Myg  = “ T«j  f(k | , e | ) (87) 

= f^,  62)  (88) 

where 

f(kn,  en)  = (X^  - \)  cos  r + (XCT  - Xg)  sin  r 

Ncl’  \l  = nondimensionai  coordinates  of  the  towpoint  with  respect  to 
keel  at  the  stem;  see  Figure  23 

Xk2-  Xe2  = coordinates  of  the  resultant  windforce 

\ = perpendicular  distance  from  the  keel  to  the  center  of  gravity 

STEADY  STATE  EQUILIBRIUM 

The  steady  state  trim  angle  r and  mean  wetted  length  Xb  are  readily  determined  from  the 
steady  state  force  and  moment  equations.  The  force  equation  is  resolved  into  a vertical  (lift) 
component  and  a horizontal  (drag)  component.  The  lift  equation  is  readily  obtained  by  setting 
the  boat  weight  W equal  to  the  sum  of  the  vertical  component  of  the  various  force  contribu- 
tions defined  earlier.  Thus  in  nondimensional  form  we  have 

W'  * - Zs'  = (Fds'  + Fgg')  cos  r - FpS'  sin  r (89) 


•Thii  it « good  iMUMpdon  for  the  model  experiment!  dbcuteed  liter. 
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where  Zs'  is  the  nondimensional  hydrodynamic  lift,  positive  down.  Similarly  the  nondimen- 
sional  tow  force  FyS'  is  set  equal  to  the  drag  or  horizontal  component  of  the  force  contribu- 
tions. 

= “ = (Fds  + F0S')  sin  r + FpS'  cos  r * Fws' 

= - Zs'  tan  r + FpS'/cos  r + Fws'  = Ds'  (90) 

where  Ds'  is  the  nondimensional  drag  of  the  boat.  The  moment  equation  is  obtained  by 
summing  the  components 

Ms'  ~ Mds'  + MBS'  + Mfs'  + Mts'  + Mws'  = 0 (91) 

where  the  component  moments  are  given  with  respect  to  the  center  of  gravity. 

If  we  substitute  Equations  (68),  (82),  and  (85)  into  Equation  (89)  we  obtain  finally  for 
the  steady  state  lift  equation 

W'  = - Zc*  = . - - sin  r cos2  t ( 1 - sin  0)  + Cn  . X sin2  r cos3  r cos  P 
a 1 + A 2 

+ x2  sin  r cos  r - X Cf  sin  t/ cos  P (92) 

Cv2 

If  we  multiply  this  equation  by  1 + X we  obtain  the  following  form  of  Equation  (92)  as  a 
cubic  equation  in  X. 

DX3  + (C  + D + E)  X2  + (B  + C - E - W')  X - W'  = 0 (93) 


where 


IT  9 

B = y sin  t cos^  t ( 1 — sin  )3) 
C=CDc  sin2  r cos3  t cos  P 
D = 0.624  sin  r cos  r/Cy2 
E = - Cf  sin  t/ cos  p 


(94) 


Substituting  Equations  (69),  <77),  (84),  (86),  (87),  (88),  and  (90)  into  Equation  (91) 


yields 


w , jt  sin  2 r ( 1 - sin  0)X 

Ms  =4  TTx 

0.624  sin 


(X  + h(r)  - Xg)  + CD  c (sin  2r)2  cos  Xg)“ 


c 2 
Lv 


-*2(fs)- 


(zs 


-XCf 


tan  r *5  - Fw<; 

cos  r cos  p 


) + Fws'  f(k2,  e2)  = 0 


(95) 
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where 


h(r)  = Xc2-|xc,  (%) 

Equations  (93)  and  (95)  must  be  solved  for  X and  r by  an  iterative  procedure  in  which 
successive  values  of  r are  assumed.  Equation  (93)  is  solved  for  X for  each  assumed  value  of 
r.  Successive  pairs  of  X and  r thus  determined  are  substituted  into  Equation  (95)  until  its 
magnitude  becomes  equal  to  zero  within  some  prespecified  amount.  In  the  present  analysis 
this  was  taken  as  0.002/Cy^. 

STATIC  STABILITY  DERIVATIVES 

As  noted  earlier  the  force  and  moment  stability  derivatives  with  respect  to  the  heave  and 
pitch  velocity  and  acceleration  perturbations  are  readily  obtained  from  Equations  (47)  and  (50). 
There  only  remains  the  task  of  obtaining  the  derivatives  with  respect  to  the  heave  and  pitch 
angle  displacement  perturbations.  These  are  readily  obtained  with  the  aid  of  the  expressions 
for  Zs'  and  M</  in  Equations  (92)  and  (95).  Since,  for  a given  deadrise  and  speed,  they  are 
functions  only  of  X(z',  r)  and  r we  have  for  the  derivatives  with  respect  to  the  heave  and 
pitch  perturbations 

az'  °ZS  ax 

0z'  " ax  az' 

ax'  _ ax*'  ax 

a7  ax  az' 

3M'  = an,’  ax 
az'  ax  az' 

az-  _ 

a o dr  + ax 

djc  + 

d6  dr  3X 

axr  _ a^s' 

a o ~ ~dr  + ax 


(97) 

(98) 

(99) 

ax 

dr 

(100) 

ax 

dr 

(101) 

ax 

dr 

(102) 

where  all  rotations  are  with  respect  to  the  center  of  gravity;  z'  = z/b  and  z equals  perturbation 
in  heave  displacement,  positive  down.  From  Equations  (90),  (92),  and  (95)  we  readily  find 


- H sin2  t cos3  r 


9ZS' 

ax 


- 2 J sin  r cos  r X + Cf  sin  r/ cos  (3 


0 03) 


9Zp'  X ^ , . 2 

1 = _ G — (cos3  t - 2 sm^  t cos  r) 

9 t 1 + X 

- HX(2  sin  t cos4  r - 3 cos2  t sin3  r) 

- J X2  cos  2 r + Cf  X cos  t/ cos  0 


9 Fee’  j 

tan  T + — Ll  — — * •— f tan  r 

9X  9X  3X  cost  2X 


9 Xs  3 Zg 


zs'  Fps'  »an  r „ »J£ ... . +JL- 

dT  ■ 3r  ,anT  + cos2T  C0ST  8t  =«2  r 


9 Xg  9 Zg 


9 Mg 
9X 


- G sin  r cos  r 


W 


(X  + h(r)  - Xg)  + y 


i_l 

+ xj 


+ H sin2  t cos2  t (X  - Xg)  + J sin  r X (X  - 2 Xg> 

cf  (<ml!  , \ f8ZS  , 9 X f ( k , , a , > 

+ STp  \-r"  V VST  “ cos  0 cos  r/  1 

^ . G -Lj  [cos  2r(X  * h(r)  - Xg)  + <0.157  - 0.00025  « j^l 
3r  1 + * L J 

+ H X ^ - Xg^  sin  2 r cos  2 r + J X2  (■—  " Xg)  cost 
\9r  cos2  r cos  0 cos  t) 

- (^'  — - - F*s')  «(k>- '•>  - Fws 


004) 

(105) 

(106) 


(107) 


(108) 


where 

g<»Ven)=  <Xkn-VsinT  -<Xen- VC0ST 
G = ir(l  - sin  P)l 2 
H = CD  c cos  $ 

J = 0.624 /Cv2 


(109) 
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Also  with  the  aid  of  Equations  (76)  and  (77)  and  Figure  20  it  is  easy  to  show  that 


3X_ i 

bz'  sir*  t 


3X 

dr 


A k - K \ 

i * + Xv)  + (0.57  + 0.001 

\ tan  r V 


0) 


tan  (3 
4 sin2  r 


(110) 

(111) 


STABILITY  EQUATIONS 


The  total  nondimensional  force  and  moment  equations  may  now  be  written  as 

'■'-T'j-V'xAY  'a  — — i d X . — i 

mu  - Ts  + Xs  + AXd  + T7  + JT  + 


m'z  = W'  + Zs'  + AZd'  + p-  + 


, . az;  az' 

z'  be 


*y  6 - Ms  +AMD  + 3z'  + 


where 


m 


r 


ATc 


m 


ATC 


JfPU2b2 


= nondimensional  mass  of  the  boat 

- nondimensional  pitch  moment  of  inertia  about  the 
center  of  gravity 

= perturbation  in  tow  force 


(112) 

(113) 

(114) 


From  the  static  equilibrium  Equations  (89)  through  (91)  we  require  that 


TS'  + Xs' 

= 0 

w'  + zs' 

= 0 

Ms' 

= 0 

(115) 

if  the  body  is  free  to  surge  then  AT^'  is  zero.  If  the  body  is  restrained  in  surge  u'  and  u’ 
are  zero,  and  the  pitch  and  heave  motions  are  uncoupled  from  surge. 

Assuming  that  the  body  is  free  to  surge,  the  stability  equations  at  the  equilibrium  flying 
conditions  are  obtained  by  putting  ATS'  equal  to  zero  and  substituting  Equations  (47),  (49), 
(50),  and  (115)  into  Equations  (112)  through  (1 14).  This  leads  to  the  following  equations 
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(116) 


where  the  primes  have  been  omitted  for  convenience. 

(X^,  - m)  u + \u  + + X*z  + Xzz  + Xg0  + Xg0  + Xg0  = 0 

Z^u  + Zuu  + (Z"  - m)z  + Z-ji  + Zzz  + Zg0  + Z g0  + Z g0  = 0 
M^u  + Muu  + M"'z  + M + Mzz  + (Mg  - ly)0  + Mg0  + Mg0  = 0 
where  the  force  derivatives  are 


Z^  = - i^(X)  sin  r cos  J n'  ds' 

X{i  = zu  tan  r 

Zu  = - 2 v?(X)  p5'  sin  r cos2  r 

Xu  = zu  tan  r 

Z~  = - i^(X)  cos2  tJ h'  ds' 

*z  = zz  tan  T 

Zz  = - 2^(X)/is'  cos3  r 

X^  = Z^  tan  r 

_ 3_£s'  ?x 
Zz " ax  az' 

X2  = Zz  tan  r 

Zg  = - .(x/mV  - s')  ds’ 

Xg  = Zg  tan  r 

Zg  = ~ 2 v?(X)Ms'  Xg  cos2  r 

Xg  = Zg  tan  r 

„ _ dZs'  dZs'  ax 
Z0  dr  + ax  dr 

Xg  = Zg  tan  r + 

and  the  moment  derivatives  are 


7 ' 

2 

cos^  r 


My=  <?(X)  sin  h'  (a'  - s')  ds' 


Mu=  2 v?(X)  sin  r cos  r 


(/"' ds'  - \ “>’) 
tJ " (a'  - s')  ds' 

= 2 ^(\)  cos2  r^/V  ds’  " 


Mj  = <f(\)  cos  r 


w _dMs'  ax 
M'~  ax  az' 


Mg=-^(X)y" Ji'(a'-s')2  ds' 


Mj=  - 2^(X)  cos  r 


Mfl- 


dMs'  dMs' 
dr  dX 


{^s  ^g2  <a’  " s ^ ds  ^ 

ax 

dv 


(117) 

(118) 


(119) 


(120) 
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The  integrations  are  from  0 to  A^.  The  evaluation  of  the  added  mass  terms  are  readily  obtained 
with  the  aid  of  Equations  (78)  through  (81 ).  The  static  derivatives  are  given  by  Equations  (97) 
through  ( 1 08).  The  expressions  for  fds’  and  Mpg'  are  given  by  Equations  (68)  and  (69). 

It  is  seen  from  Equation  (119)  that  for  small  trim  angles  the  derivatives  of  the  X-equation 
are  considerably  smaller  than  those  of  the  Z-equation.  Under  these  conditions  it  is  reasonable 
to  expect  that  the  influence  of  the  surge  degree  of  freedom  on  the  porpoising  stability  would 
be  small  enough  so  that  the  X-equation  could  be  omitted.  In  any  case,  if  the  boat  is  being 
towed  at  constant  speed,  then,  as  previously  noted,  values  of  u and  u in  Equations  (117)  and 
(118)  are  zero.  This  effectively  uncouples  the  surge  degree  of  freedom,  and  the  stability  may 
be  determined  from  the  pitch  and  heave  equations  alone.  In  the  present  investigation,  all  of 
the  data  available  for  checking  the  theory  were  obtained  by  using  planing  boat  models  towed 
at  constant  speed.  For  these  reasons  numerical  calculations  have  been  made  using  only  the 
pitch  and  heave  stability  equations. 
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APPENDIX  B 


ESTIMATES  OF  EFFECT  OF  WINDAGE  AND  CHINE  RADIUS 
ON  STABILITY  DERIVATIVES 

A high  metal  framework  was  attached  to  the  models  of  Day  and  Haag.*  Its  purpose  was 
to  support  a scribe  for  recording  the  motion  after  the  boat  started  to  porpoise.  This  introduced 
a small  aerodynamic  drag  and  moment  on  the  model.  The  following  estimated  values  for  this 
drag  were  used  in  the  computations 

Fws'  = 0.0032 
Xe2  = 1.3 
Xk2  = 2.6 

Shuford21  found  that  a small  radius  on  the  chine  will  account  for  a reduction  in  lift  of 
5 to  10  percent  - corresponding  to  a 1/64-  and  1/1 6-inch  radius,  respectively,  on  a 4-inch-beam, 
flat-bottom,  planing  surface.  On  the  basis  of  chine  radius  measurements  obtained  from  a 
typical  wood  model,  it  was  estimated  that  the  models  used  by  Day  and  Haag  had  about  a 5 
percent  loss  in  dynamic  lift.  The  theoretical  calculations  were  therefore  modified  to  take  this 
into  account.  This  amounted  to  reducing  the  magnitude  of  P and  Q in  Equations  (78)  through 
(81)  and  G and  H in  Equations  (103)  through  (108)  by  5 percent. 
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Figure  I - Coordinate  System 
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EQUILIBRIUM  TRIM  ANGLE, 


Figure  3 — Variation  of  Stability  Root  with  Trim  Angle  for  Various  Values  of 
Speed  Coefficient,  Deadrise  Angle  of  0 Degree,  Load  Coefficient  of  0.72 
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Figure  4 - Comparison  of  Theoretical  and  Measured  Porpoising 
Boundaries  for  Deadrise  of  0 Degree 


39 


2 


3 


4 


5 


, DEGREES 


Figure  6 - Comparison  of  Theoretical  and  Measured  Porpoising 
Boundaries  for  Deadrise  of  20.5  Degrees 
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Figure  7 - Variation  of  Critical  Trim  Angle  wi*h  Load  Speed  Factor 
for  Various  Loadings.  Deadrise  Angle  of  0 Degrees 


42 


CRITICAL  TRIM  ANGLE,  r„,  DEGREES 


16 


14 


12 


8 


6 


4 


2 

0.1  0.2  0.3  0.4 

LOAD  SPEED  FACTOR,  s/CLb/2 


Figure  9 - Variation  of  Critical  Trim  Angle  with  Load  Speed  Factor 
for  Various  Loadings.  Deadrist  Angle  of  20.5  Degrees 
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CRITICAL  TRIM  ANGLE,  r . DEGREES 


Figure  9 - Variation  of  Critical  Trim  Angle  with  Load  Speed  Factor 
for  Various  Loadings,  Deadrist  Angle  of  20,5  Degrees 
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CRITICAL  LONGITUDINAL  CENTER  OF  GRAVITY  TO  BEAM  RATIO,  X 


Figure  10  - Comparison  of  Measured  with  Computed  Position 
of  Center  of  Gravity  at  Porpoising 
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CRITICAL  MEAN  WETTED  LENGTH  TO  BEAM  RATIO.  \ 
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Figure  1 1 - Variation  of  Mean  Welled  Length  to  Beam  Ratio  at 
Porpoising,  with  Speed,  Deadrise  Angle,  and  Loading 
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Figure  12  - Variation  of  Critical  Trim  Angle  with  Load  Speed 
Factor  and  Radius  of  Gyration-Beam  Ratio 
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Figure  13  - Variation  of  Critical  Trim  Angle  with  Load  Speed  Factor  and  Load  Coefficient 
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Figure  14  - Effect  of  Radius  of  Gyration  and  Speed  on  Location  of  Center 
of  Gravity  at  Porpoising  (Tow  point  at  center  of  gravity) 


CRITICAL  LONGITUDINAL  CENTER  OF  GRAVITY  TO  BEAM  RATIO,  X 


Figure  15  - Effect  of  Load  Coefficient  and  Speed  on  Location  of  Center 
of  Gravity  at  Porpoising  (Tow  point  at  center  of  gravity) 
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Figure  16  - tffcct  of  Rad:,,s  of  Gymtion  and  Speed  on  Mean  Welled 
Length  at  Porpoising  (Tow  point  at  center  of  gravity) 
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52 


') 


VOLUMETRIC  FROUDE  NUMBER,  U>/gV1/3 


Figure  18  - Comparison  of  Approximate  Estimates  of  Porpoising  Boundaries 
of  Reference  8 with  DTMB  Series  62  Model  Data 
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ire  19  — Comparisons  of  Approximate  Estimates  of  Pm<vnsing  Boundaries 
of  Reference  I 1 with  DTMB  Series  62  and  EMB  Series  5G  Model  Data 


Figure  20  — Geometric  Relationships  in  Steady  State  Planing 
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Figure  21  — Correction  Factors  for  Three-Dimensional  Flow 


Figure  22  — Variation  of  Flow  Coefficient  with  Deadrise  Angle 
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Figure  23  - Coordinate*',  of  Tow  Point  and  Center  of  Wind  Force 
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